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Abstract 

Improving performance in motor skill acquisition is proposed to be supported by tuning of neural networks. To address this 
issue we investigated changes of phase-amplitude cross-frequency coupling (paCFC) in neuronal networks during motor 
performance improvement. We recorded intracranially from subdural electrodes (electrocorticogram; ECoG) from 6 patients 
who learned 3 distinct motor tasks requiring coordination of finger movements with an external cue (serial response task, 
auditory motor coordination task, go/no-go). Performance improved in all subjects and all tasks during the first block and 
plateaued in subsequent blocks. Performance improvement was paralled by increasing neural changes in the trial-to-trial 
paCFC between theta {9; 4-8 Hz) phase and high gamma (HG; 80-180 Hz) amplitude. Electrodes showing this covariation 
pattern (Pearson's r ranging up to .45) were located contralateral to the limb performing the task and were observed 
predominantly in motor brain regions. We observed stable paCFC when task performance asymptoted. Our results indicate 
that motor performance improvement is accompanied by adjustments in the dynamics and topology of neuronal network 
interactions in the 9 and HG range. The location of the involved electrodes suggests that oscillatory dynamics in motor 
cortices support performance improvement with practice. 
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Introduction 

Phase-amplitude cross-frequency coupling (paCFC) of oscilla- 
tions in different frequency bands has been proposed as an 
effective mechanism to form functional networks that recruit local 
neuronal populations across a global spatial scale [1-4]. Phase- 
amplitude CFC between HG (80-150 Hz) amplitude to 9 (4- 
8 Hz) phase was first described by [2] and later confirmed by other 
authors in rats [5,6] and humans [7] . During paCFC amplitudes 
of higher frequency oscillations, reflecting local cortical processing, 
are modulated by the phase of low frequency oscillations [8-12]. 
This mechanism has been proposed to engage and coordinate 
local processing modules across spatially distributed brain areas 
supporting cognition and motor performance [4,9, 13- 17]. Further 
support for this proposal comes from recent clinical studies linking 
altered paCFC to debilitating psychiatric and motor disorders 
[18-21]. Moreover, paCFC is prominent during language and 
motor tasks [2,4] and the frequency of the slower phase coupling 
oscillation is task dependent [12]. However, beyond clinical studies 
evidence for a functional role of paCFC in the process of 
organizing human cognition and behavior is limited predomi- 
nandy to the memory domain (see [22] for a review). Axmacher 
and colleagues [23] reported that inter-individual differences in 



working memory performance correlated with differences in 
paCFC precision, supporting the functional relevance of CFC 
for memory processing. Tort and colleagues [6] examined the 
dynamic modification of functional relations between performance 
and CFC in rat hippocampus and found coupling strength 
between 9 and gamma (y: 25-100 Hz) correlated with maze 
learning. 

A stronger link between paCFC and behavior in humans would 
be supported by a correlation between paCFC and trial-by-trial 
variations in performance. To address this, we examined the 
relation between paCFC and motor performance improvement. 
We recorded the electrocorticogram (ECoG) in human patients 
(N = 6; mean age = 20.5, std =5.5; 2 female) undergoing epilepsy 
diagnosis while they learned skilled motor behaviors. To assess the 
link between paCFC and behavior we compared changes in 
paCFC to changes in performance over an extended time scale 
during motor skill acquisition, and correlated performance and 
paCFC at the single trial level. We show that paCFC in 
intracranial subdural recordings between 9 (4—8 Hz) and HG 
(80-180 Hz) in the human cortex tracks level of motor perfor- 
mance across different motor tasks. 
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Results 



experimental session 



Phase amplitude cross frequency coupling 

We investigated potential links between paCFC and motor 
performance in six subjects each performing one of three repetitive 
motor tasks described next. The three different behavioral tasks 
(Figure 1) were a serial reaction time task (SRT 2 subjects), a 
go/no-go task (GNG 3 subjects), and an auditory motor 
coordination task (AMCT 1 subject). All three tasks required the 
coordination of finger movements with an external stimulus. We 
assessed motor performance as reaction time in the SRT and 
GNG tasks and as the temporal deviation from the target time 
point in the AMC task. The cognitive requirements for 
performance improvement are different in all three tasks: learning 
the motor sequence in the SRT, learning the stimulus-response 
association in the GNG and improving movement timing in the 
AMCT. However, the motor component is performance improve- 
ment with practice. The dynamics of the different performance 
measures were assessed in a group statistic, by comparing the 
average behavioral outcomes between fixed trial bins (see Figure 2). 
We recorded the ECoG while subjects performed two blocks of 
one of each task with the hand contralateral to the electrode grid. 
The ECoG-time series were filtered in the 9 -band (4-8 Hz) and in 
the HG-band (80-180 Hz) yielding two separate filtered signals 
(see Methods). We calculated the analytic amplitude of the HG- 
band time series by taking the absolute value of the Hilbert 
transform of the filtered time series. The analytic amplitude is a 
new time series representing the amplitude envelope of the HG- 
oscillations at any moment in time. We performed the analysis on 
the 500 ms interval immediately following the stimulus onset. This 
interval includes the preparation of the responses indicated by the 
stimulus and includes approximately three 0-cycles. 

We first asked whether the amplitude envelope of the local HG 
oscillations is phase coupled to the local 0-band oscillations. 
Figure 3A shows the time course of sine waves fitted to the single 
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Figure 2. Here we depict the separation of the whole 
experimental session into trial bins. The experimental session in 
each patient consisted of 2 blocks separated by a short break. In each 
block we defined two trial bins each containing 30 trials (blue). We 
compared the PLV across the four trial bins to assess the evolution of 
connectivity length of and HG activity during motor performance 
improvement. 

doi:10.1371/journal.pone.0089576.g002 

trial variations of HG analytic amplitude pooled across all 
electrodes in one subject. As predicted, HG analytic amplitude 
varied systematically over the 0-cycle. Figure 3B shows single 
subject sine wave fits to the HG analytic amplitude averaged over 
trials and electrodes. Each fit was significant (p< 0.001) and the 
HG amplitude variations were consistent over subjects with only a 
slight deviation of subject 1 (see Figure 3B). The frequencies of 
the fitted sine waves are in the 8 band (.95 Hz, SE: .02 Hz) and 
the phase angle is .6 rad (SE: .23 rad, see Methods for an 
explanation of the sine wave parameters - frequency and phase 
angle). The maximum of the HG-analytic amplitude centers 
around the trough in the 0-cycle (mean = 2.56 rad, std = .56 rad, 
skewness = —.16). 
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Figure 1. Paradigms employed (details described in Methods). A) Serial reaction time task: The numbers on the screen indicate the finger to 
be used for the key press. B) Go/no-go: Green indicates a go and red indicates a no-go trial. C) Auditory motor coordination: Subjects were instructed 
to press a key in the middle of the interval between two consecutive tones. The interval length was either one second or two seconds and was held fixed 
for one minute. Each subject carried out two blocks (see Methods). 
doi:1 0.1 371 /journal.pone.0089576.g001 
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Figure 3. The amplitude of the HG oscillations is phase coupled to the 0 -band (4-8 Hz) oscillations in all subjects across paradigms. 

A) Time courses of sine wave functions fitted to the single trial amplitude envelopes of the HG oscillations of one subjects collapsed over electrodes. 

B) Sine wave functions fitted to the trial-averaged HG oscillation amplitudes envelopes of each subject. Each solid line represents the fit for one 
subject. Each dot represents the individual trial average of the HG oscillation in one of 20 intervals equally spaced over a 0 cycle. The black dashed 
line shows the averaged sine waves across subjects. The vertical blue dashed line denotes the averaged phase angle the HG amplitude peaks across 
subjects. The maximum of the 0 cycle is at phase 0 and the minimum at +n. 

doi:1 0.1 371 /journal.pone.0089576.g003 



Covariation of paCFC with performance improvement 

To investigate if cross frequency coupling (paCFC) covaries with 
motor performance variations, we first calculated the trough to 
peak ratio (^-trough to HG peak ratio - TPR; see Methods) over 
all electrodes on the grids as a metric for paCFC and related it to 
behavioral performance. Figure 4A shows the development of 
TPR and motor performance over the time course of the two 
experimental blocks each subject completed. Both TPR and 
performance increased during the experiment, as indicated by the 
fitted exponential functions. A statistical test confirmed this 
finding. In this test, we first compared average motor performance 
in the initial 30 trials of the first block with performance in the final 
30 trials and found a significant improvement (Wilcoxon rank sum 
test across all subjects: p < .05, Figure 4B, See Figure 2 for the 
structure of the experimental session and Table 1 for mean 
performance measures for each trial bin). However, performance 
plateaued in the second block as indicated by no significant 
difference (p = 0.18). The difference between the first and the 
second block is indicated by a significant block-by-trial-bin 
interaction in a two way ANOVA across subjects (F(l,20) = 
1 1.28; p = .003, df error = (Nsubj -lfNtrialbin). The next question 
was whether TPR exhibits the same behavior (Figure 4C). In 
concordance with behavioral performance we found, that the TPR 
increased between the first and the last 30 trials of the first block 
(Wilcoxon rank sum test across all subjects: p < .05) but did not 
change between the first and the last 30 trials of the second block 
(p = 0.3). A significant block-by-trial-bin interaction in a two way 
ANOVA (F(l,20) =5.95; p = .03) confirmed that TPR changed 
during the first block and plateaued during the second block. This 
suggests that, on average, paCFC covaries with motor perfor- 
mance with paCFC and motor performance increasing early in 
the first experimental block and both plateauing in the second 
block. 

Support for a functional relation between paCFC and motor 
performance would be provided by a trial-by-trial TPR with 
performance correlation. In order (i) to test for this correlation and 



(ii) to disentangle cortical regions showing varying paCFC with 
performance, we pooled the data in six anterior and frontal 
regions of interest (ROIs); the anterior and posterior medial frontal 
gyrus, the anterior and the posterior inferior frontal gyrus, and the 
superior and inferior sensorimotor cortex (see Figure 5) for all five 
subjects with a square 8x8 (N = 4)/16xl6 (N = 1) grid implanta- 
tion (see Figure S2). In each ROI we pooled the TPR values across 
electrodes and determined the p-values of the trial-by-trial 
correlation with performance of each ROI (Figure 5 for details 
see Methods and Figure SI). Significant correlation of TPR with 
motor performance (corrected for multiple comparisons) was 
observed in pre-/motor cortex and in anterior and posterior 
inferior frontal sulcus. We predicted two sources of variability in 
single trial correlation between TPR and motor paCFC: one that 
is performance improvement related and varies systematically over 
time and another one that is not related to performance 
improvement and varies randomly from trial to trial. The first 
analysis supported performance improvement related trial-by-trial 
correlations reflecting the co-evolution of coordination between 
brain networks and improvements of motor performance. We then 
calculated in the same ROIs the partial correlation of TPR with 
performance. This analysis factored out the fraction of correlation 
between TPR and performance which can be attributed to 
random trial-by-trial covariations and is performance improve- 
ment unrelated. This performance improvement unrelated corre- 
lation of TPR with motor performance was observed in 
sensorimotor cortex and in premotor cortex, in the posterior 
middle temporal sulcus (corrected for multiple comparisons) and 
overlaps with the performance improvement related correlation. 

Discrimination of performance improvement from 
Random Performance Fluctuations 

To disentangle these two potential and functionally distinct 
causes of paCFC-performance covariations and disentangle the 
spatially wide ROIs we performed two different correlation 
analyses separately for each recording electrode. First, we 
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Figure 4. Covariation of average paCFC with performance over the time course of the experiment. A) The development of TPR and 
motor performance during the time course of the experiment. Results are collapsed across all six subjects/three experiments. Data for the first block 
and second block are shown in the left and second half of the plot. Each point represents the average in one of 1 00 time bins. Exponential functions 
fitted to the data z-scored over both blocks indicate a similar time course for performance and TPR. B) Subject averaged motor performance during 
the first and last sets of 30 trials in the first (early learning) and the second (late learning) experimental block. C) Subject averaged TPR. Data was z- 
scored within blocks and TPR was averaged over all electrodes. 
doi:1 0.1 371 /journal.pone.0089576.g004 



calculated the partial correlation of TPR with performance. 
Second, we calculated the standard Pearson correlation between 
TPR and motor performance. Pearson correlation captured the 
performance improvement related plus the performance improve- 
ment unrelated correlation. Combined with partial correlation this 
was used to distinguish between the two effects. Electrodes that 
capture performance improvement related TPR with performance 
correlations should show a partial correlation different from zero 
and a Pearson correlation different from the partial correlation (see 
Methods). Importantly, we reasoned that if we observed a 
significant Pearson correlation in an electrode that significantly 
changes if we discount time related correlations (in partial 
correlation), then the TPR - performance correlation in this 
electrode is partly due to performance improvement related TPR 
- performance correlations. Figure 6A shows the electrodes where 
random trial-by-trial fluctuations of TPR correlated with motor 
performance (significant partial correlation - uncorrected for 
multiple comparisons). Clusters of electrodes showing high 



Table 1. Behavioral data. 





Paradigm 


trial bin 








Patient 


1 


2 


3 


4 


SRT 


SRT01 


917 (271) 


767 (227) 


748 (219) 


707 (161) 


SRT02 


1472 (300) 


1017 (284) 


859 (207) 


966 (232) 


AMCT 


AMCT01 


1 1 7 (66) 


98 (65) 


147 (197) 


142 (129) 


Go/No-Go 


GNG01 


343 (46) 


323 (117) 


331 (109) 


302 (128) 


GNG02 


426 (185) 


301 (42) 


260 (28) 


301 (66) 


GNG03 


433 (168) 


379 (178) 


286 (40) 


393 (89) 



For SRT (serial reaction time task) and GNG (Go/No-Go) task reaction time is 
shown (standard deviation) in msec. For AMCT (auditory-motor coordination 
task) the absolute deviation from precision is shown also in msec. Each trialbin 
encompasses 30 trials. 
doi:1 0.1 371 /journal.pone.0089576.t001 



correlations are located in sensorimotor cortex, in premotor 
cortex, in lateral prefrontal cortex and in ventral anterior temporal 
cortex. Figure 6B shows the distribution of electrodes with 
performance improvement related trial-by-trial correlations be- 
tween TPR and motor performance. Clusters of performance 
improvement related electrodes were apparent in premotor cortex, 
in lateral prefrontal cortex and in ventral anterior temporal cortex. 
Importandy, the variation of TPR with performance improvement 
was not a result of a shift of the HG amplitude peak relative to the 
8 trough and hence the coupling phase remained stable during 
performance improvement (see Appendix SI). 

Discussion 

Phase-amplitude cross-frequency coupling has been proposed to 
support interaction within functional networks [4] . Here we show 
that fluctuations of 6 and HG paCFC are tightly linked to motor 
performance improvement at the single-trial level and show robust 
performance improvement clusters over pre-motor and motor 
cortices. 

Performance improvement and theta and high gamma 
activity 

During motor performance improvement 9 and HG activity 
show a clear development of coupling that asymptotes in strength 
as motor behavior performance improvement plateaus. paCFC is 
highly dynamic and task-specific [4,24] and it has been proposed 
that paCFC enables adaptive behavior [2]. Here we report data 
from three different motor behavior experiments showing that 
paCFC dynamics reflect adaptive behavior supporting a relation- 
ship between paCFC and motor performance improvement on a 
trial-to-trial level. Notably, despite the differences in tasks similar 
cortical regions associated with performance improvement or 
random trial-to-trial performance are identified by paCFC 
evolution. The dynamic nature of paCFC and the task-specific 
coupling patterns have been shown in a variety of studies with task 
dependent differences in coupling frequencies and coupling phase 
[12,24—27]. Here, we add an important paCFC characteristic. We 
show that even though the preferred phase as indicated by 
coupling phase stability does not change the activity pattern of 
both frequencies varies with behavioral changes. 
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Figure 5. Depiction of the results from the ROI-analysis. A) ROIs with significant performance improvement unrelated TPR/performance 
correlations. B) ROIs with significant performance improvement related TPR/performance correlations. ROIs with significant correlations (Bonferroni 
correct for six comparisons) are marked with an asterisk. The 6 ROIs are the anterior and posterior medial frontal gyrus, the anterior and the posterior 
inferior frontal gyrus, and the superior and posterior sensorimotor cortex. The blue margin shows the grid coverage across all subjects with a square 
grid implanted. 

doi:1 0.1 371 /journal.pone.0089576.g005 



Biological mechanism 

Oscillatory dynamics are proposed [1] to be inherent to the 
interplay of brain regions for cognitive control in memory and 
learning [21]. For example 9 activity observed in hippocampal and 
neocortical regions varies as a function of the state of the subjects 
[28]. The neocortex exhibited more prominent 9 activity during 
wakefulness compared to REM sleep. Performance improvement 
during practice can be achieved by distributed 9 networks - which 
are up-regulated during wakefulness - by integrating or coordi- 
nating local activity. Here, the concept of information integration 
means that 6 oscillatory activity accumulates and integrates the 
results of local processing as reflected in HG activity in the 
premotor/motor region. HG activity, either an indication of 
spiking activity or very fast network oscillations [25,29], may be 
involved in planning and initiation of motor responses [30] . This 
frequency possibly reflects the activation of cortico-subcortical 
networks involved in the feedback control of discrete movements 
[31]. Taken together we speculate that information on planning of 
motor responses is integrated into memory by paCFC which 
results in performance improvement during the process of 
performance improvement. 

Conclusion 

We identified cross-frequency coupling in the human cortex 
which is associated with motor performance variability per se. In 
this network a smaller area is integrated whose oscillatory 
dynamics reflect the progress in performance improvement. This 
learning related network suggests the establishment of a memory 
trace which is accumulated during practice and which is 
represented in a mutually adapted level of activity of 9 and HG 
activity [16,32-34]. In this respect paCFC provides a mechanism 
subserving motor memory formation [2]. 

Materials and Methods 

Patients 

Six epilepsy patients undergoing pre-surgical monitoring with 
subdural electrodes participated in the experiments after providing 



their written informed consent. Experimental and clinical record- 
ings were taken in parallel. Recordings took place at the University 
of California San Francisco (UCSF), CA, USA (4 Patients), Johns 
Hopkins University, Baltimore, USA (1 Patient) and the Epilepsy 
Center Bethel (ECB), Bielefeld (1 Patient), Germany and were 
approved by the local ethics committees ("Committee for the 
Protection of Human Subjects at UC Berkeley", "Johns Hopkins 
Medicine Institutional Review Board" and "Ethical Committee of 
the University of Magdeburg"). 

Experimental Paradigms 

We carried out three different motor tasks (serial reaction, go/ 
no-go, auditory-motor coordination) with six different patients 
(Fig. 1). Each patient participated in one of the tasks. All 
paradigms required coordination of key presses on a computer 
keyboard to an external stimulus. Patients performed the task 
sitting upright in their bed using the hand contralateral to the grid. 

Serial Reaction Task. The serial reaction task (SRT) 
consisted of a series of visually cued finger taps. The subjects 
had their fingers placed on different keys of a laptop keyboard 
(right hand: space bar, j, k, and ; - left hand: space bar, f, d, 
and a). Trials started with one of the numbers 1, 2, 3, or 5 
appearing on a laptop-screen cueing the movement of thumb, 
index finger, middle finger, or little finger, respectively. Numbers 
were presented until a key was pressed but maximally for 2 seconds. 
In each subject the four numbers were presented in a fixed 
sequence (six items long) or random order depending of the block 
number but only fixed blocks were used. Each block took 
approximately 10 minutes. Two patients participated in this task 
(SRT01 - 02). They were instructed to press keys as fast and 
accurate as possible. One block took approximately six minutes. 

Auditory-Motor Coordination Task. The second motor- 
paradigm was an auditory-motor coordination task (AMCT). One 
patient participated and was instructed to respond as accurately as 
possible halfway between successive auditory clicks presented at a 
constant rate. Seven click sequences, each 60 s long, were 
presented in a block. The inter-click-interval in a sequence was 
either 500 ms, 1000 ms or 2000 ms and the participant was 
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Figure 6. Electrodes with significant trial-by-trial correlations of TPR with performance. The significance threshold was determined in a 
permutation procedure (see Methods) A) Learning unrelated correlations of TPR with performance. B) Learning related correlations of TPR with 
performance. Darker colors indicate stronger correlations. See Methods for calculations on separating performance and learning related effects. The 
blue shape in the first and second row show the outline of all superimposed square grids. The black shapes in the third row denotes the grid 
locations for the participant in the AMCT. Spatial distortions result from the projection onto the cortex (for details see Figure S2). 
doi:1 0.1 371 /journal.pone.0089576.g006 



informed that the interclick-interval changes. The seven click 
sequences with the differing interclick-interval were presented 
randomly within each block. Only sequences with a interclick- 
interval of 1000 and 2000 ms entered the analysis. The clicks were 



placed in front of the patient at 1 m distance. One block lasted 
7 min. 

Go/No-Go task. Three patients participated in the go/no-go 
(GNG) task (GNG01 - 03). In each trial the patients were 



presented with speakers plugged into the laptop's sound card and presented with either a green or red square of 100 ms duration 
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and 900 - 1 200 ms stimulus onset asynchrony. The subjects were 
instructed to respond as quickly as possible to green squares by 
pressing a key on a laptop keyboard but to withhold responses 
when a red square was presented (in block 1-3: 20% of the trials; 
in block 4: 50% of the trials). The participants were familiarized 
with the task in an initial short practice session. Only correct Go- 
trials of the first two blocks entered the analysis. Each block lasted 
approximately eight minutes. 

Data recording 

At UCSF the electrocorticogram (ECoG) was recorded either 
from 64 platin-iridium-electrode grids arranged in an 8x8 array 
with 10 mm center-to-center spacing (FTT01, FTT02, GNG01, 
GNG02 GNG03) or from a 256 electrode grid (both Ad-Tech 
Medical Instrument Corporation, Racine, Wisconsin) arranged in 
a 16x16 array with 4 mm center-to-center spacing (GNG02). 
Exposed electrode diameter was 2.3 mm in the 64 electrodes grid 
and 1 .8 mm in the 256 electrodes grid. The electrode signals were 
recorded with a 256 channel preamplifier (PZ2-256, Tucker-Davis 
Technologies (TDT), Inc) with the electrode furthest from the 
motor cortex used as a reference. The data from the pre-amplifier 
were sampled at 3051.7 Hz on a digital signal processor (RZ2 4 
DSP, Tucker-Davis Technologies (TDT), Inc) with 16-bit resolu- 
tion and stored to hard disk. Trigger signals indicating button 
presses and stimulus onsets were sent from the stimulus laptop via 
a USB-1208FS DAQ. (Measurement Computing, Norton, MA) 
plus a photodiode attached to the screen and recorded on the DSP 
synchronized to the brain data. Trigger timing was additionally 
recorded on the stimulus laptop by querying the computers 
performance counter using the Psychophysics Toolbox (www. 
psychtoolbox.org). In Bielefeld (AMCT) the ECoG signal was 
recorded at 1000 Hz sampling frequency ( 1 6 Bit resolution) with a 
Nihon Kohden system (Tokyo, Japan) equipped with auxiliary 
analogue channels for synchronous recording of the trigger signals 
and the output from the sound card. Here 5 stripes were 
implanted each equipped with two parallel rows of 5 electrodes 
each (see Figure 6). 

Data analysis 

We used Madab 2008a (Mathworks, Natick, USA) for aU offline 
data processing. We first preprocessed the recorded brain data and 
then we derived measures quantifying adaptation of oscillatory 
neural dynamics during motor skill learning. All filtering was done 
using IIR filters (Butterworth filter of order 4). Preprocessing 
served to remove non-physiological artifacts from the recorded 
data and to prepare them for further analysis. First we excluded 
channels exhibiting ictal activity or excessive noise from further 
analysis. In the remaining good channels we then excluded 
time intervals containing artifactual signal distortions such as steps 
of pulses by visual inspection. Finally, we re-referenced the 
remaining electrode time-series by subtracting the common 
average reference 



1 

xcAn(t)= -y~)x c (t) 

c= 1 



(1) 



calculated over the n good channels c from each channel time 
series. The resulting time series were then used to characterize 
brain dynamics over the time course of motor behavior 
performance improvement in terms of the TPR: 

For each trial starting at stimulus onset we calculated the TPR 
to quantify the evolution of phase-amplitude cross-frequency 
interactions of cortical oscillations during motor skill learning. 



Therefore, we band-pass filtered each electrode's time series at 
two frequency bands, in the 0-range (4—8 Hz) and in the HG (80- 
180 Hz) range since coupling was task relevant between these 
frequencies across a variety of experimental tasks [2] . We detected 
^-troughs, the local minima, in the 0-range filtered time series in 
the interval between 0 to 500 ms after stimulus onset (Figure 7). 
We obtained the HG analytic amplitude Anoit) by Hilbert- 
transforming the HG filtered time series. For each detected 6- 
trough we then estimated the depth of the through Dg and the 
simultaneous HG amplitude as the average of the ^-filtered and 
the Aug time series over an interval of 83 ms (half 0 oscillation) 
centered on the trough. Note that multiple 8 troughs fit into the 
500 ms analysis leading to multiple estimates per trial. We 
averaged the individual estimates Dg{t) and Ajj(}(t) to obtain 
one measure for 9 trough depth Dg and one for HG amplitude 
Ahg for each trial j. From these values we calculated TPR for 
each trial j as: 



TPR f = log 



(2) 



HGJ 



Taking the log of the ratio makes the distribution of TPRs 
symmetric. Note that the TPR includes both stimulus-locked and 
non-stimulus-locked brain activity. It summarizes the global cross 
frequency interaction on the grids. 

Phase-Amplitude Coupling 

paCFC was tested by splitting the 6 oscillations of the 500 ms 
analysis window into 20 equally spaced phase bins ranging from -n 
to 7t (18°or 0.314 rad) in each subject and each electrode. In each 
phase bin we averaged the amplitude envelope of the local HG. A 
cosine wave function 
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Figure 7. Calculation of the trough to peak ratio (TPR). We 

quantified paCFC as the ratio between 0 trough (local minima of the 0 
time series - red vertical lines) and HG amplitude at the corresponding 
0 trough. Around each detected trough we spanned a window (half 0 
cycle - gray bars) in which 0 activity (black bold line) and HG amplitude 
(green line) was averaged. 
doi:10.1371/journal.pone.0089576.g007 
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y = a * cos (co * HG + 0) (3) 

with a representing the amplitude, ft) representing the frequency 
and <j> representing the phase angle was fitted to the resulting 20 
HG amplitude (y) values, co close to 1 indicates that HG amplitude 
variation is accounted for the 8 cycle. 

ROI analysis 

In each patient we grouped electrodes according to the same 
anatomical landmarks in 6 regions of interest (see Figure SI): The 
anterior (in sum 34 electrodes across subjects) and posterior (44 
electrodes) medial frontal gyrus, the anterior (49 electrodes) and 
the posterior inferior frontal gyrus (8 1 electrodes), and the superior 
(48 electrodes) and inferior (46 electrodes) sensorimotor cortex. We 
averaged the TPR-values within each ROI across electrodes. In 
each patient we determined the p-value for Pearson's correlation 
coefficient r and the partial correlation coefficient rho (p) between 
the averaged TPR values and behavioral performance across 
trials. Each ROI in which the mean p-value across subjects fell 

.05 

below the p-value corrected for multiple comparisons (p C orr = —7-) 

6 

was considered statistically significant. 

Separating performance from learning effects 

We separated performance from learning effects by applying a 
permutation test statistic. The reasoning for applying a permuta- 
tion test was two-fold. First, we sought to correct the p-values for 
each electrode due to the many individual correlation tests 
applied. We tested this against a distribution which did not rely on 
the same temporal interval (500 ms following the stimulus 
presentation) for which the correlation coefficient was calculated. 
Second, we wanted to identify electrodes in which Pearsons 
correlation coefficient r was significandy higher than the partial 
correlation coefficient p. This means that we looked for electrodes 
with a significant difference between r and p. Since the 
significance can only be determined in relation to a distribution 
we estimated this distribution from our data. Hence, the null 
hypothesis to be rejected was that the difference between 
electrodes r and p was derived from a random distribution. The 
recorded time series were filtered in the 6 (4—8 Hz) and in the HG 
(80-180 Hz) frequency. Subsequendy, we calculated the HG 
envelope of the HG time series in each electrode and each trial by 
taking the absolute value of the Hilbert transform of the filtered 
time series. The analytic amplitude is a new time series 
representing the amplitude envelope of the HG-oscillations at 
any moment in time. We then determined 20 time windows 
around the stimulus onset each with a width of 500 ms and 
400 ms overlap. In order to conduct the TPR permutation test 
statistic we calculated Dg and Aug (see above) around each 9 
trough in the time window in each electrode and trial which yields 
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